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Abstract 

The main physics programme of the International Linear Cohider (ILC) requires a 
measurement of the beam energy at the interaction point with an accuracy of 10~^ 
or better. To achieve this goal a magnetic spectrometer using high resolution beam 
position monitors (BPMs) has been proposed. This paper reports on the cavity 
BPM system that was deployed to test this proposal. We demonstrate sub-micron 
resolution and micron level stability over 20 hours for aim long BPM triplet. We 
find micron-level stability over 1 hour for 3 BPM stations distributed over a 30 m 
long baseline. The understanding of the behaviour and response of the BPMs gained 
from this work has allowed full spectrometer tests to be carried out. 
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1 Introduction 



The physics potential of the International Linear Collider (ILC) depends 
greatly on precise energy measurements of the electron and positron beams at 
the interaction point (IP). Two types of analysis are particularly sensitive to 
the collision energy: threshold cross-section measurement and reconstruction 
of particle resonances [Tj. The required accuracy for the mass measurements 
dictate that the fractional error on the determination of the incoming beam 
energy must be better than 10~^. To measure the energy to this level and 
to minimise the disruption of the beam, a magnetic spectrometer has been 
proposed. 

When passing through a magnetic field, a particle with an electric charge q is 
deflected by an angle 9 which is inversely proportional to its energy E\ 



where c is the speed of light, B is the magnetic field and d£ is the path segment 
along which the particle travels. The initial ILC spectrometer proposal [2] 
suggested using a single bend and, with careful characterisation of the bending 
magnets and an accurate measurement of the bend angle 9, the energy of the 
beam could then be reconstructed. 

The performance of a similar spectrometer has already been demonstrated 
during the second phase of the Large Electron Positron (LEP2) collider at 
CERN P|. The magnetic spectrometer installed at LEP2 achieved an accu- 
racy of 1.7 X 10"'^ by measuring the change in bend angle as the beam passed 
through a single steel dipole magnet. To obtain the high accuracy required, 
a relative energy measurement was made with the spectrometer calibrated at 
the resonance using the resonant depolarisation method, thus removing the 
need for an absolute angle measurement. In order to further improve the accu- 
racy of the measurement, an average was taken over many bunches and many 
revolutions of the beam around the accelerator. This allowed the requirements 
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Fig. 1. Schematic of the current basehne design for the ILC spectrometer, 
on the BPM resolution to be relatively loose. 

As the ILC is a linear machine, resonant depolarisation is not an option for 
the beam energy measurement and in order to limit low energy operation, the 
ILC spectrometer will have to provide an absolute energy measurement. In 
addition, a bunch-by-bunch energy measurement is highly desirable to remove 
some of the dependence on the stability of the machine during physics data 
taking. These constraints, in addition to the required resolution of the energy 
measurement (10~^), force the performance of the combined BPM and elec- 
tronics to be considerably higher than the system used at LEP2. In order to 
achieve this level of performance, high resolution cavity BPMs will need to be 
used. However, as cavity BPMs are sensitive to the beam slope as well as the 
beam offset, we are focusing our studies on a different approach to the energy 
measurement problem from the single bend method initially proposed. Two 
identical magnets with opposite fields allow the beam energy to be measured 
as a function of the induced horizontal displacement x given by 



where d is the total bend length and the deflection angle 9 is small (fig. [T]). 
The inclination of the beam trajectory through the BPMs is then minimised 
and, by mounting the mid-chicane BPMs on high precision movers, they could 
be translated to the beam location in the case of large induced offsets. The 
precision and accuracy of the measured offset x contributes directly to the 
uncertainty of the final energy measurement. The current baseline spectrome- 
ter design ^ implements a 5 mm deflection. An offset measurement to better 
than 500 nm is therefore required to achieve the necessary energy resolution. 
However, a BPM system with better performance would allow a smaller deflec- 
tion and, therefore, reduce the emittance growth due to synchrotron radiation. 
In addition, the system has to be stable to the level of 500 nm over multiple 
hours to avoid extensive recalibration and consequent loss of luminosity. 

The operation and stability of a long baseline BPM system in the proposed 
spectrometer is of interest in several other fields. For example, the beam po- 
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Fig. 2. Schematic of the SLAC beam hnes that provide a high energy electron beam 
from the 2-mile SLAC linac, through the Beam Switch Yard to the A-line leading 
to End Station A. 

sition in the ILC linacs must be measured with a resolution of 1 fim for the 
orbit correction and emittance preservation. High resolution BPMs are also 
required throughout the beam deUvery system (BDS). Current test facilities 
and modern linac based light sources are also demanding significantly better 
performance from their BPM systems. Most notably, the ATF2 [5] facility at 
KEK, designed to test the final focus optics for the ILC, not only requires 
resolutions of ~ 100 nm for the extraction line BPMs but good stability and 
ease of use of the entire monitoring system. 

In the framework of the T-474 test beam experiment [6] in End Station A 
(ESA) at the Stanford Linear Accelerator Center (SLAC), we installed several 
BPM stations in the 40 m available drift space. The aim of three running 
periods in 2006 was to commission and optimise these BPMs and study their 
resolution and stability. The principal results of these runs are discussed below. 



2 ESA beam and hardware configuration 

2.1 Beam delivery to End Station A 

SLAC's accelerator and beam transport system for delivering high energy 
electron beams are depicted in Figure [2l A high intensity electron beam is 
produced by a thermionic gun, bunched and pre-accelerated in the first sec- 
tions of the linac to 1.19 GeV. The bunches are then kicked by a pulsed 
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Table 1 

Beam parameters at ESA and as proposed for the ILC [7]. 



Parameter 


SLAG ESA 


ILC-500 


Repetition rate 


10 Hz 


5 Hz 


Energy 


28.5 GeV 


250 GeV 


Train length 


up to 400 ns 


1 ms 


Micro bunch spacing 


20 - 400 ns 


337 ns 


Bunches per train 


1 or 2 


2820 


Bunch charge 


1.6 X 10^° 


2.0 X 10^° 


Bunch length 


500 /im 


300 /um 


Energy spread 


0.15% 


0.1% 



magnet into the Linac-to-Ring transfer line and then transported to the elec- 
tron North Damping Ring (DR), where they are stored for 8 ms to reduce the 
beam emittance. The Ring-to-Linac transfer line transports the bunches back 
from the DR and a pulsed magnet kicks them into the linac at Sector 2. The 
beam is subsequently accelerated to 28.5 GeV. At the end of the linac it is 
then transported from the Beam Switch Yard (BSY) through a 24.5° bending 
section into the A-line on its way to End Station A. ESA test beams operate 
with single bunches at 10 Hz parasitically to PEP-II operation. The ESA is 
currently the highest energy test beam facility available, with its other beam 
parameters, listed in Table [H similar to those for the ILC. 

The A-line leading from the BSY to ESA is about 300 m long (see Figure [3]). 
Two BSY dipoles, Bl and B2, just upstream of the D-10 beam dump at the 
start of the A-line bend the beam through an initial angle of 0.5°, and are 
followed by a string of 12 dipoles which deflect the beam further by 24°. The 
quadrupole lattice in the A-line consists of a doublet, QIO and Qll, at the start 
of the A-line which brings the beam to a horizontal waist at the dispersion 
matching quadrupoles, Q19 and Q20. These are located at the high dispersion 
point (see Figure H]) in the middle of the A-line to correct first- and second- 
order horizontal dispersion in ESA. Four additional quadrupoles (Q27, Q28, 
Q30, Q38) near the end of the A-line control the transversal beam size in ESA. 
The beta functions and horizontal dispersion through the A-line and ESA are 
shown in Figure HI Vertical and horizontal corrector dipole pairs, A28/A29 
and A32/A33 are used to set the beam trajectory in ESA. They are used in 
the SLAC Control Program (SCP) feedback to stabilise the beam position at 
BPM 31 at the end of the A-line and BPM 1 in ESA. SL-10 is a high power 
momentum slit located at the A-line bend mid-point where the dispersion 
is at maximum. The maximum momentum acceptance is 2%. A SCP energy 
feedback system uses a stripline BPM, BPM17, at a high dispersion point near 
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Fig. 3. The beam line configuration for the A-line. The beam goes from left to right, 
top to bottom. 



SL-10 to stabilise the beam energy. 



A-line beam diagnostics include RF cavity BPMs, charge-sensitive toroids, 
a synchrotron light monitor, retractable profile screens, and high frequency 
diodes. The synchrotron light monitor is positioned upstream of the mid-bend 
region, and is used to monitor beam energy spread and provide a diagnostic 
for minimising it. High frequency diodes are installed in ESA for monitoring 
and tuning the bunch length [S]. 



2.2 ESA beam line and experimental equipment 



The configuration of beam line components, beam diagnostics and experi- 
mental equipment at the end of the A-line and in ESA is shown in Figure O 
Two protection collimators are located in ESA; 3C1 (19 mm aperture ra- 
dius) is at the entrance to ESA in front of BPM 1 and 3C2 (8 mm aperture 
radius) is located in front of BPM 3. There are two beam profile monitors, 
one upstream of the T480 collimator wakefield experiment (PR2) and one 
just beyond the east wall of ESA (PR4). These are aluminium oxide screens 
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Fig. 4. The beta functions {Px, Py) and horizontal dispersion (Dx) in the A-line and 
ESA. The origin is the swivel point in the BSY. The maximum dispersion is at the 
SL-19 momentum slits in the Aline, as well as the horizontal waist. 
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Fig. 5. Schematic of the principal beam line components at the end of the A-line 
and in ESA. 



that can be inserted into the beampipe. Two wire scanners, WSl and WS2, 
are approximately 20 m apart and are used to perform transverse beam size 
measurements. 75 /im tungsten wires are scanned across the beam, creating 
low energy electrons that are detected by photomultiplier tubes at 90°. Several 
beam loss monitors along the beam line in ESA are interlocked to the machine 
and radiation protection systems. 

Our experimental equipment includes a BPM doublet (BPMs 1,2) and two 
BPM triplets (BPMs 3-5 and 9-11) and an interferometer system monitoring 
horizontal motion of the BPM 3-5 triplet. These systems are described in 
more detail below. Among other experimental setups, a collimator wakefield 
experiment is located downstream of the BPM 1 and 2 doublet and uses the 
BPM system for measuring wakefield kick angles [S]. 

We measured the beam jitter at various BPM stations along the beam line. 
The typical values for the data discussed in this paper are detailed in Table El 
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Table 2 

The typical beam jitter as recorded at the BPM stations for the data considered 
in this paper. The error indicates the variation seen over a one hour run. 



BPMs 


Jitter X ( fim) 


Jitter y ( ^m) 


1-2 


110 ± 17 


14± 11 


3-5 


157 ± 15 


40± 19 


9-11 


173 ± 16 


66± 22 



2.3 Cavity BPMs 



A microwave or Radio Frequency (RF) cavity BPM is a discontinuity in the 
beam line which, through the excitation of different oscillating electromagnetic 
field configurations (cavity modes), can be used to measure the position of the 
transiting bunch. It can be shown that the amplitude of the excited dipole 
field, providing the bunch is not far from the centre of the BPM, is not only 
proportional to the charge, but also related to the beam offset x from the 
centre, the slope 6 of the beam trajectory and the tilt a of the bunch with 
respect to the cavity [TOpTj : 

Vx oc xe~'"* sin ut 

Va oc —ae~^^ cos tut (3) 
Vg oc 6e^^^ cos cut 



where the V^s are contributions to the induced voltage, and uj and P are the 
frequency and the decay constant of the dipole mode. The combined cavity 
output is given by: 

V{t) = Vx + V^ + Ve (4) 
= e"^* [A^x sin(cjt) -|- {Ag9 + A^a) cos(cjt)] , 



Therefore the position component is always in quadrature with the combined 
slope and tilt signal and consequently, the angle and position components can 
be separated providing the beam arrival time is known. For that purpose and 
also as a bunch charge reference we use additional smaller cavities operating 
with the monopole mode at the same frequency as the dipole mode in the 
position cavities. Their signals follow the same processing path as the BPM 
signals in order to preserve the phase relation. 

We used two types of cavity BPMs in our studies: rectangular and cylindrical, 
all of which have a similar dipole mode frequency. The loaded quality factor 
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Table 3 

Parameters of the various BPMs along the ESA beam hne. The distances in brackets 
given for the location are measured relative to the entrance of the experimental hall. 



XJX J.VJ. 1J.IU.J.J.1IL/ ^ J. 




T^pson^irit frpn 

J. l/tyOWl.J.Clji.J.lJ J-J. VjVJ • 


TiOJ^npn (J T 


A ■nprfnrp 






(MHz) 




(mm) 


31, 32 


A-Line 


2856 


~ 3000 


51 


1, 2 


ESA (~ 10 m) 


2856 


~ 3000 


51(1) and 38(2) 


3, 4, 5 


ESA (~ 25 m) 


2859 


~ 500 


36 


9, 10, 11 


ESA (~ 40 m) 


2856 


~ 3000 


20 



Ql, which defines the signal decay constant T via 



is also similar for all the cavities except for 3, 4 and 5. The general cavity 
parameters are summarised in Table [31 

Seven rectangular cavity BPMs are available: two in the A-line (BPMs 31 
and 32) and five in the End Station (BPMs 1, 2, 9-11). BPMs 9, 10 and 11 
were originally designed for use in the main SLAC linac whereas BPMs 31, 
32, 1 and 2 were built for the A-Line. Each BPM consists of three separate 
cavities (see Figure [6^): one cylindrical reference cavity and two rectangular 
cavities delivering x and y position dependent signals. The rectangular shape 
of the position cavities helps reduce the effects of coupling between the x and 
y orientations. We tuned all the rectangular BPMs (including their monopole 
cavities) to a nominal frequency of 2.856 GHz using external tuners. These 
BPMs were initially installed for the SLAC E158 experiment |12] . During this 
installation, BPM 1 suffered significant mechanical damage and as a conse- 
quence, both the X and y cavity signals had strong coupling with other modes 
resulting in worse resolution and stability. 

The central BPM station (BPMs 3-5) consists of three cylindrical cavities 
designed for use in the cryogenic regions of the ILC linac (see Figure [GJd) ^T3\ . 
In these BPMs, a single cavity is used to provide both x and y signals. A 
combination of slots and waveguide couplers provides very good suppression 
of the unwanted monopole modes. These cavities have a lower Q-value and 
therefore a shorter decay time of the dipole mode signal to provide bunch- 
to-bunch position measurements in the ILC without the need to exclude the 
contamination from the previous bunches. We used the reference cavity of 
BPM 9 to provide the phase and charge information for BPMs 3-5; it was 
therefore tuned to 2.859 GHz. 
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Fig. 6. a) An example of the rectangular cavities used for BPMs 31, 32, 1, 2, 9, 10 
and 11 showing x, y and reference channels and b) the central BPM of the prototype 
ILC linac BPMs triplet. 

2.4 DAQ and signal processing electronics 

The signal processing electronics that convert the 2.9 GHz signals coupled out 
of the cavities to low frequency signals prior to digitisation consist of a single 
stage down-mixing circuit (see Figure [Tj). The amplitude of the input signal can 
be adjusted manually with variable attenuators and phase adjustment is also 
possible. The signal is filtered at the front end with a narrow 20 MHz band- 
pass filter centred at 2856 MHz to suppress the noise and unwanted modes. No 
amplification is applied at this stage in order to use the full dynamic range of 
the mixers. The filtered signal is then mixed with a local oscillator (LO) signal 
operating at 2939 MHz to produce a low frequency waveform at ~ 83 MHz. An 
image reject mixer scheme is used to reduce the impact of the unwanted noise 
and cavity modes around the image frequency. The down-mixed signal is then 
passed through 20 dB of amplification and an additional 20 MHz bandpass 
filter. 

The electronics for BPMs 3-5 and 9-11 are located close to the beam line, 
about 20 m from the BPM stations in the experimental hall. The electronics 
for BPMs 31, 32, 1 and 2 are in a different area ("Counting House"), with an 
additional ~ 30 m of cabling required to transfer the high frequency cavity 
signals to the electronics racks. Also, the signals from BPMs 31, 1 and 2 
are split to provide input for the SLAC Control Program (SCP) running the 
position feedback, adding 6 dB attenuation to these channels. 

Digitisation is carried out by four VME Struck SlS330x digitisers, two 12 bit 
(BPMs 31, 32, 1 and 2) and two 14 bit (BPMs 3-5 and 9-11), with ±2.5 V 
input range and 50 Vt input impedance. An external 119 MHz clock signal 
derived from the linac RF system is used for all the modules. Consequently, 
the digitiser under-samples the signals, producing a waveform aliased to ~ 
36 MHz. Under-sampling is employed in order to widen the gap between the 
signal peak and image sideband, thus simplifying the hardware filtering at the 
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Fig. 7. Schematic view of the signal processing electronics. The components within 
the dashed lines are contained in a single mixer module 

front end. The digitisers are triggered using a pulse synchronous to the bunch 
arrival time, which is generated by the linac control system. 

The data from the SIS digitisers, as well as that from additional gated Analog- 
to-Digital Converters (ADCs), VME Smart Analog Monitor (VSAM) and 
Computer Automated Machine and Control (CAMAC) modules are captured 
by a Windows XP PC running Labview DAQ programme . In addition to the 
waveform data, some environmental data are also available such as informa- 
tion on the linac status from the SCP, temperature data from thermocouples 
placed on the BPMs and electronics, interferometer data (see Section 12. 6p . 
stray magnetic field data measured with the flux gate probe and high voltage 
data for the wire scanners and other experiments. Data are stored at varying 
rates depending on the source: the waveform data are recorded on a per event 
basis (10 Hz) with each event corresponding to one bunch being delivered 
from the linac, but data provided by the SCP and other additional modules 
are recorded at a significantly slower rate (~ 0.1 Hz). 



2. 5 Digital signal processing 



The voltage at the front-end of the analogue electronics given by Equation <^ 
can be re-written as: 

V{t) = Ae~^'?>m{ut + (l)) (6) 



where A = ^ [A^xY + {A06 + A^aY is the amplitude of the waveform and 
= tan-i is the phase. 

For simplicity, we can assume that the analogue electronics are ideal and only 
convert the frequency of the waveforms to the first intermediate frequency 
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~ 83 MHz so that for both the position and reference cavities we have: 



V{t) = ae"^* sm{uJiFit + (f)) (7) 

Vreiit) = Orefe"^* Sm{uJiFlt + 0ref) (8) 



where we neglect the differences of the resonant frequencies and decay con- 
stants between the cavities. Note that the decay constant of the signals is 
changed by the band pass filters. 

The waveforms described by Equations ([7]) and ([HD are what appears at the 
front-end of the digitisers. Due to under-sampling, the frequency of the wave- 
form experiences another conversion to //i?2 = fs — fiFi, where fs is the 
sampling frequency of the digitisers (119 MHz as mentioned above). The am- 
plitude of the resulting waveform remains the same. 

In order to extract the amplitude and phase information necessary to recover 
the position, this waveform is down-converted again in software by multiplying 
by a complex LO signal at the same frequency as the waveform. This process 
is called digital down-conversion (DDC). It results in a signal describing the 
envelope of the initial waveform: 



DDC 



(t) = ae-^* sin (w/^a^ + ^)e-''^'^^^ (9) 



2 



The mixing process also produces the up-converted component with twice 
the initial frequency. In order to remove this component and further suppress 
out-of-band noise, a digital filter is applied. We chose a Gaussian filter as it 
is easy to implement and has a good phase response within the pass band. 
Filtering is done in the time domain. As we are interested in the amplitude of 
the waveform at only one sample point, we reduce the computation time by 
restricting the filtering to one summation out of the whole convolution: 



=+k 



^DDC(^O) = E ■ F{tn - to) (10) 



where V^ddc is the filtered voltage response after down-conversion at the time 
to, k defines the window of the filter, outside which the contribution is ne- 
glected (we set the limit to 0.1% of the maximum filter response), and the 
filter function is given by: 

, , V2^Af /-t2(27rA/)2\ ^ ^ 
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Fig. 8. A down-converted waveform as recorded at the digitiser as well as the am- 
plitude and phase response of the DDC algorithm for a whole waveform. 

where A/ is the bandwidth of the filter and fs is the sampling frequency. 

Applying the filter across all samples we can compute the envelope of the am- 
plitude and the phase for the whole waveform (see Figure [8]). A smooth enve- 
lope with no oscillations and a flat phase response in the region corresponding 
to large beam power are an indication of a good frequency determination and 
low contamination of the waveform with any unwanted signals. The phase vari- 
ation seen prior to this region is a consequence of the broad-band components 
of the initial beam excitation. 

The filtered DDC response depends on the amplitude a and phase of the 
original digitised waveform: 



Vi^cito) cx ae'^' (12) 

^DDC,ref(^0) OC a^efC^^^^ (13) 



where (j)' includes any phase delays occurring. As mentioned above, the position 
can be extracted with the help of a reference signal. Normalising the complex 
amplitude from the position channel with that from the reference channel, we 
exclude the charge dependence of the position signal: 

^DDcfa) ^ ^gi(<A'-<^;ef) = —e'^'l' (14) 

DC, re f (to) Ctref Oref 



where the phase offset A0 between the position and the reference signals at 
the time to does not change unless some changes occur in the processing chain. 

In order to recover the position information, we need to calibrate the BPMs 



13 



Table 4 

The relative misalignment of the BPMs with respect to the mean straight line 
orbit where the error quoted is the error on the mean. The alignment of BPM 4 was 
optimised continually through the run as it was mounted on the mover system. 



BPM 


Offset X (nm) 


Offset y {^m) 


1 


-133 ±3 


168 ±0 


2 


-90 ±3 


-94 ±0 


3 


-128 ±4 


-193 ± 1 


5 


-134 ±5 


10 ± 1 


9 


59 ±6 


-60 ±2 


10 


125 ±6 


29 ± 2 


11 


302 ±6 


140 ±2 



applying a known offset either to the beam or the BPM. From this we measure 
the scale S (which is the sensitivity of the term ^ to the beam offset) and 
the phase ^/q (IQ - in-phase/quadrature rotation) of all the points with a 
positive offset. The position is then given by: 



X 



S Re 



«ref 



(15) 



Thus, with an accurate calibration the beam position can be found from the 
amplitude and phase of a chosen sample of the waveform. The calibration 
coefficients 5* and $/q will only stay accurate for a limited period of time due 
to gain and phase drifts in the processing electronics and cables, frequency 
changes of the cavities and other environmental factors. We therefore took 
calibration data several times during the running period. 



In order to optimise the performance of the BPMs and electronics, both the 
attenuation and alignment offsets of the BPMs needed to be minimised. Ini- 
tially, with the attenuation set close to maximum to avoid saturation of the 
waveforms, we measured the relative misalignment of the BPMs from the nom- 
inal beam trajectory. Using these data we realigned the BPMs, enabling us to 
reduce the attenuation. We then repeated this process several times to opti- 
mise the sensitivity and dynamic range. The final offset values of the BPMs 
relative to the mean straight line orbit are shown in Table |H 
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Fig. 9. The interferometer system in ESA. a) A schematic of the laser path in 
the interferometer relative to BPMs 3-5 and the beam line, b) The interferometer 
elements mounted on the optical table with BPMs 3-5 on the right. The laser path 
is indicated schematically. 

2. 6 Mover system and interferometer 

The middle BPM of the central triplet (BPM 4) is mounted on a dual axis 
mover system that allows travel in x and y directions. The stages each have 
a 2 mm pitch lead screw driven by a 200 steps per revolution stepper motor. 
This gives a resolution of 10 /xm per step for the BPM moves. The motion 
system additionally features an LVDT position readout with an accuracy of 
~ 6 /xm. The maximum allowed travel range of the both stages is limited to 
about ±6 mm for radiation and machine protection reasons. For calibration, 
a travel distance of ±500 yum is used, thus giving a contribution of ~ 1% to 
the uncertainty on the calibration scale. 

To provide information on the exact x position of BPM 4 and hence allow 
corrections to the calibration and mechanical position of the assembly, we used 
a set of three single-pass linear interferometers [H] manufactured by Zygo. 
Each interferometer measures the relative horizontal displacement of the BPM 
with respect to interferometer heads on an aluminium table (see Figure |9]). 
A single heterodyne helium-neon laser supplies all three interferometers with 
polarised light. The interferometer heads are polarised beam splitters that 
provide the reference and target arms of each of the interferometers with one of 
the two laser frequencies. Photo-detectors on a Zygo 4004 VME measurement 
board extract the phase and amplitude of the optical beat frequency present 
in the recombined light path, allowing a measurement of the relative velocity 
of the BPMs. The measurement board provides a single bit position resolution 
of approximately 0.3 nm. 
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Table 5 

Total motion and non-rigid body (NRB) motion measured over 30 minutes in the 
X direction of BPMs 3-5. The errors quoted are the spread of the RMS. 



BPM 


RMS (Total) 


RMS (NRB) 


3 


170 ± 4 nm 


94 ± 1 nm 


4 


680 ± 37 nm 


620 ± 10 nm 


5 


130 ± 3 nm 


72 ± 1 nm 



Using the interferometer we measured the vibrational motion of BPMs 3-5. 
This not only includes the rigid motion of the entire system, but also the non- 
rigid body motion (i.e. the deviation of each of the individual BPM positions 
from a straight line) which contributes directly to the measured resolution of 
the BPMs. The amount of motion recorded at each BPM as well as the non- 
rigid body motion as extracted from the interferometer data using an SVD 
method (described in Section 14. ip are listed in Table [5l As BPM 4 is mounted 
on a non-rigid support, this shows significantly higher vibrational motion than 
the other two BPMs. 



3 Commissioning of the BPMs 

Section m will be devoted to the system's performance during an 18 hour long 
continuous data taking period which took place at the end of the 2006 running, 
while in this section we will first report on our studies of the methodology of 
dealing with BPM signals. More specifically, we will discuss the extraction of 
the calibration coefficients described in Section [23| namely the BPM frequency 
a;, the phase difference between position and reference signals ^jq and the 
scale S, as well as the optimisation of the DDC algorithm in order to obtain 
the best resolution and stability of the BPMs. 

For the discussion here and the rest of the analysis results we excluded BPM 
pulses that were saturating the digitisers due to inaccuracies in determining 
the extrapolation factor. 

3. 1 BPM frequency and sampling point 

The frequency of the digitised BPM signal needs to be measured to a consid- 
erable degree of accuracy in order to ensure phase stability when performing 
the DDC (see Section 12.51) . To determine the optimum frequency for each 
cavity, we fitted the waveforms to an exponentially decaying sine wave using 
MINUIT [15]. To avoid low amplitude effects, this fit was only performed on 
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Table 6 

The frequencies that were used throughout the 18 hour run. The BPMs with errors 
of zero had a neghgible systematic error. The frequency for BPM llx is for the first 
8 hours only. 



BPM 


fx (MHz) 


fy (MHz) 


fQ (MHz) 


1 


36.10 ±0.06 


39.47 ±0.03 


36.70 ± 0.02 


2 


36.66 ±0.01 


36.83 ±0.01 


as ql 


3 


38.46 ± 0.03 


37.73 ± 0.08 


39.16 ±0.01 


4 


39.38 ± 0.07 


38.85 ± 0.04 


as q3 


5 


38.94 ±0.15 


39.10 ±0.04 


as q3 


9 


36.53 ±0.01 


36.65 ±0.01 


as qlO 


10 


36.65 ±0.01 


37.29 ± 0.00 


37.05 ±0.00 


11 


36.46 ± 0.00 


36.84 ±0.01 


as qlO 



events where the beam was driven to large offsets (~ 0.5 mm, which corre- 
sponds to waveforms of approximately half the dynamic range of the digitiser). 
We combined the fits from all such events over an 18 hour period and took the 
median in order to get an accurate estimation of the frequency for each BPM 
throughout the run. The uncertainties associated with these measurements 
are entirely systematic, due to either physical drift and/or systematic errors 
associated with the fitting method. We therefore defined the uncertainty on 
the frequency as the RMS of its distribution. The frequencies determined in 
this way are listed in Table [6l 

The variation in frequencies seen over the 18 hour period was caused by two 
different effects. BPMs 3-5 had significantly larger variation than the majority 
of the remaining BPMs due to systematic effects in fitting waveforms with a 
short decay time. The number of samples that could be included in the fit was 
significantly lower than for the other BPMs, thus leading to a larger error. 
The variation seen in BPMs 9-11 was consistent with being caused by thermal 
expansion of the cavities. The frequency of a cavity BPM is roughly linearly 
dependent on the size of the cavity. Assuming a thermal expansion coefficient 
for copper of 1.7 x 10^^ per °C, a change in temperature of ~ 0.1°C leads 
to an expansion of 1.7 x 10~^, equivalent to ~ 5 kHz at S-band. Typical 
frequency variation over the 18 hour period is shown in Figure [TOk. The BPM 
frequency showed a correlation with temperature (see Figure [TOb) which is 
in reasonable agreement with this prediction. Despite having a similar decay 
time, the variation for BPM 1 was significantly higher than BPMs 9-11 and 
2. This was found to be due to interference between different cavity modes as 
a result of the mechanical damage to the cavity. 

Variations in frequency compromise the calibration and produce an appar- 
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Fig. 10. a) The median frequency determined from large offset (~ 0.5 mm) events 
over the 18 hour running period for BPM lOy. b) The fitted frequency versus BPM 
temperature. 

ent drift and a reduction in resolution if not re-calibrated. Our studies have 
shown however that induced drifts and resolution deterioration due to fre- 
quency changes are small compared to those introduced by other effects (see 
Section m). 

During the 18 hour period, there was a significant change in the frequency of 
BPM llx from 36.46 MHz to 36.55 MHz. This change occurred during a short 
access to the experimental hall, so it seems most likely that this was caused 
by a mechanical disturbance of the cavity or the associated tuner. During the 
remainder of the 18 hour running, the altered frequency was applied where 
appropriate. 

After we determined the frequency, we could then find the sampling point {to 
in Equation [TOj) . This point is defined as the location in the down-converted 
waveform where we sample the phase and amplitude for the BPM considered. 
For obvious reasons, the Gaussian filter window is centred around this sam- 
ple. There is a number of considerations when choosing the optimal sampling 
point for the analysis including the maximisation of the signal-to-noise ra- 
tio, ensuring the sampling point is away from any monopole contamination 
at the beginning of the waveform and also the phase being constant around 
the selected sample. We decided to pick a point 10% down from the peak of 
the filtered waveform (see Figure [8]), as a trade-off between the criteria men- 
tioned above. We found the optimum sampling point for each BPM analysing 
the data containing a large range of the offsets and taking the mean of the 
resultant sampling point distribution. 

Variation of the time difference between the digitiser trigger and the beam 
arrival time would cause several systematic errors. Firstly, any difference in 
frequency between the position and reference cavities would make the relative 
phase sensitive to trigger time variation, thus changing the IQ phase. Though 
we had attempted to minimise this effect by tuning both the reference and 
position cavities to the same frequency, differences were still present. Secondly, 
the trigger time variation would change the relative position of the fixed sam- 
pling point with respect to the waveform. This would cause variations in both 
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Table 7 

The filter bandwidth used in the DDC analysis for each set of BPMs. 



BPM 


Bandwidth ( MHz) 


1-2 


1.3 


3-5 


3.8 


9-11 


1.3 


Q cavities 


2.8 



the recorded amplitude as well as the phase (if the phase was not stable around 
the chosen sampling point) that would produce apparent changes in scale and 
IQ phase. As the digitiser was triggered by signals synchronised with the linac 
timing system, its variation was negligibly small. We determined such varia- 
tion to be less than a fraction of the digitiser sample by measuring the timing 
jitter of a reference cavity waveform rectified with a crystal detector. 

The only remaining parameter to optimise was the bandwidth of the Gaussian 
filter used in the DDC algorithm (A/ in Equation [TT]) . To determine the 
optimal value for the filter bandwidth, we calibrated the BPMs and calculated 
the resolution (as defined in Section 14.11) for a range of bandwidth values. For 
each BPM station, we used the bandwidth that gave the best resolution for the 
remainder of the tests (see Table [7]). Note that the optimal filter bandwidth 
scales with the inverse of the cavity fill time, as expected. 

3.2 Calibration 

To measure the remaining calibration coefficients, S and $/q in Equation [T5| 
we induced a known change in the recorded beam position either through 
movement of the beam or the BPM. The magnitude of these moves is restricted 
by the dynamic range of the BPM on the one hand and the beam jitter and 
short time scale drift on the other hand. We chose typical ranges of ±500 /im. 

3.2.1 BPM calibration with mechanical movers 

The most accurate method of calibration available was by means of the mover 
system on which BPM 4 is mounted. By moving this central BPM of the triplet 
and using data from the two surrounding BPMs, it is possible to remove the 
beam jitter and drift from the calibration scan in the central BPM to first 
order. These three BPMs are also instrumented in the x direction with an 
interferometer which allowed accurate calibration of the mover system in this 
axis. As the accuracy of the interferometer is ~ 5 nm, the large moves of 
±500 /zm could be established with a negligible error. 
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The method used to cahbrate BPM 4 with such data was developed by the 
NanoBPM collaboration [IT]. The beam orbit through the three BPMs (3-5) 
can be considered a straight line in both xz and yz planes. The position in 
the central BPM can then be predicted from a linear combination of the outer 
BPM coordinates: 

Xk = Sk + '^{aiXi + PiUi + pix'i + aiHi) (16) 



where Xi and i/i are the positions recorded in BPM i with the primed coordi- 
nates indicating the beam slope values and a, /?, p, a and 6 are constants that 
encode the relative rotations, offsets, scales and IQ phase differences between 
the BPMs. As these should not change over short time periods, the same con- 
stants can be used to predict the position for many events. Consequently, this 
relationship can be rewritten in terms of a matrix multiplication: 



^ Xkl ^ 



Xk2 



\XkjJ 



1 Vii x\^ ■ ■ 

1 Xi2 yi2 X'i2 V'i2 ' ' 

\ 1 Xi^ x\^ ■■■) 



fSk\ 

ai 

A 
pi 

V ■■ J 



or Xk = M • Kk 



(17) 



where the i sub-index denotes BPMs 3 and 5 and the j sub-index denotes 
event number. Thus, to find the constants a, (3, p, a and 6, the matrix of 
positions above must be inverted: 

Kk = ■ Xk (18) 



To perform the inversion, we employed the method of Singular Value Decom- 
position (SVD). This involved decomposing the (possibly singular) M x N 
matrix into the product of an M x iV column-orthogonal matrix, an x 
diagonal matrix with positive or zero elements (the singular values) and the 
transpose of an A^ x A^ orthogonal matrix which could then be trivially in- 
verted. The benefit of this method is that, as no exact solution to Equation [TSl 
is possible due to its over-constrained nature, the least squares minimisation 
of the equation is found instead. Thus, this method can be used to minimise 
the BPM residuals by finding the best set of coefficients that can predict the 
position in the central BPM from the surrounding BPMs. 

This method can also provide predictions for the in-phase (/) and quadrature 
{Q) components of the signal. First, the positions and slopes of the outer BPMs 
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must be replaced with their measured Is and Qs in Equation [181 Then the 
deconvolution and inversion must be carried out twice: once with the position 
vector (xk) replaced by a vector of the measured in-phase values of the central 
BPM and the second time with the position vector being replaced by a vector 
of measured quadrature values. 

Assuming the BPM response is linear, the prediction from the outer BPMs 
to the central BPM is independent of any beam motion that is seen by all 
BPMs. This includes both beam jitter and beam drift. However, if the central 
BPM is shifted from its position relative to the outer BPMs, the coefficients 
would still predict the position at the central BPMs original location. Thus, 
by comparing the measured and predicted Is and Qs, the consequence of this 
move can be seen. 

To perform a BPM calibration using this method, we moved the central BPM 
from —0.5 mm to +0.5 mm in 5 steps monitoring the x-moves with the in- 
terferometer. We then determined the optimum set of coefficients to predict 
the Is and Qs at the central BPM using data from the middle step in the 
calibration (i.e. at zero relative offset). From these coefficients we predicted 
the Is and Qs for the rest of the calibration. The mean values for I and Q 
at each mover position in the calibration were found by fitting the / and Q 
distributions with a Gaussian function. In the IQ plane these points form a 
straight line (see Figure [TTi) . As the source of the movement was known to be 
produced by a position change, this line therefore describes the position axis 
and consequently, its angle of inclination is equal to the IQ phase ($/q). The 
resulting predicted position in BPM 4 can then be easily calibrated against 
the known mover positions, yielding the scale S. In order to avoid saturation 
at the extreme mover positions, we only used the central three steps of the 
five step mover scan to determine the calibration constants. 



3.2.2 Corrector scan 

To calibrate the remaining BPMs, we varied the integrated field in the correc- 
tors A32 and A33 (see Figure [3]) from —0.01 kGm to +0.01 kGm in five steps. 
This induced a change of between ±0.5 mm and ±1 mm at the downstream 
ESA BPMs. We took 60-70 events per step and, for each of these, computed 
the Is and Qs. Using the change in / and Q for the various steps, we calculated 
the IQ phase ($/q) and scale {S) as for the mover calibration. Again, only 
the central three steps were used to avoid saturation. 

A significant drawback of this calibration method is that the beam jitter and 
drift can not be removed by using outer BPMs to predict the position of a 
central BPM. Though this has little impact on the IQ phase determination, 
the scale determination is very sensitive to beam drifts, which alter the as- 



21 



0) 

in 

(0 

Q. 

C 



10000 



10500 



11000 11500 
Event Number 



2 -0.05 

T3 
(0 

o 



10000 



10500 



11000 11500 
Event Number 




£ -0.05 - 
■O I 

(0 : 



In Phase Comp. /A.U. 




Predicted Position /mm 



Fig. 11. A typical mover calibration. The change of a) the In-Phase (/) and b) the 
Quadrature (Q) components against the event number and c) against each other, 
d) The predicted position against measured position after the calibration. 
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Fig. 12. A corrector scan calibration for BPM lOx. a) The In-Phase (/) and b) the 
Quadrature {Q) components against event number, and c) against each other, d) 
The predicted position against measured position after the calibration. 

sumed relationship between the true beam position and the magnetic field 
in the corrector magnets. In addition, the accuracy of the predicted beam 
positions was not known as the magnetic field in the corrector magnets was 
not continuously monitored. The actual changes in beam position could not 
therefore be determined with sufficient accuracy and this directly contributed 
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to the uncertainty on the scale. 



To estimate the improvement in scale determination when using the movers 
as well as estimate the non-linearities involved, we fitted the predicted versus 
measured position distribution (see Figure fTTH and Figure [T2li) with a straight 
line and found the maximum deviation from the fit. This gave a residual of 
~ 5% for the corrector scan and ~ 0.5% for the mover scan. 

In order to minimise the uncertainty on the scale from drifts during the scan, 
the more accurate scale for BPM 4 as found from a mover scan was used to 
correct the BPM scales found from the corrector scan: 



where Si is the scale for BPM i computed using a corrector scan, S'4,mover is 
the scale for BPM 4 found using a mover scan and SI is the corrected scale. 
The typical correction factor was of order 15%. 



4 Experimental results 

Having described the methodology of our experiment in the previous sections, 
we will here report on resolution and stability measurements over short and 
long periods of operation for the separate BPM stations as well as the complete 
eight-BPM system over a baseline of 32 m. 

4-1 BPM resolution 

To measure the resolution of a BPM, we again used the SVD method described 
in Section [3. 2. II in order to minimise the effect of any calibration inaccuracy on 
the measurement. The width of the distribution of residuals however contains 
contributions from all BPMs considered. Considering Equation [16] and making 
the assumption that the contributions from both the slope and the offset in 
the orthogonal axis are small we write the width as: 



Si — Si X 



s 



'4, mover 



(19) 




(20) 



where Rk is the resolution of BPM k. 
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BPM 10x Predicted /mm BPM 10x Residual liim 



Fig. 13. a) The predicted position versus measured position for BPM Wx and b) 
the residual distribution for the same data. 

Table 8 

The measured and predicted resolutions of the BPMs in the ESA beam line. 



BPM 


Measured x 


Predicted x 


Measured y 


Predicted y 




(/im) 


(fim) 


(fim) 


(fim) 


1, 2 


1.1 ±0.1 


1.3 ±0.7 


2.2 ±0.3 


1.4 ±0.8 


3, 4, 5 


0.53 ±0.05 


0.2 ±0.1 


0.46 ± 0.02 


0.3 ±0.2 


9, 10, 11 


0.19 ±0.01 


0.1 ±0.1 


0.17 ±0.01 


0.2 ±0.1 


Linked 


0.82 ±0.14 




1.19 ±0.35 





With the further assumption that the predicted position in the residual calcu- 
lation is defined purely by geometry, the resolutions of all BPMs are the same 
and the BPMs are equi distantly spaced. Equation [20] reduces to Rk = ^j2/3ak 

for a triplet and = ^Jlj2ak for a doublet. This correction was applied in 
the calculation of the individual BPM resolution inside a BPM station. More 
BPMs can be included, up to the entire set, and the resolution of the recon- 
structed orbit can be determined. As some of the assumptions above are not 
valid in this case, we quote the width of the distribution of residuals "as is" 
without any corrections and refer to it as the "linked" system resolution. In the 
energy spectrometer, the precision of the orbit reconstruction in the middle 
of the chicane directly contributes to the energy measurement. We therefore 
calculated the residual for BPM 3 in the middle of the baseline versus the rest 
of the system. The resolutions as recorded during a typical run are shown in 
Table [8] and an example of a residual distribution is shown in Figure [131 

The main factors contributing to the resolution are electronic noise, digitiser 
errors and mechanical motion. To estimate the combined RF electronics and 
digitiser contributions, we ran the DDC algorithm over waveforms that did not 
contain beam data (i.e. that contained only noise). Dividing the mean of these 
results by the mean reference amplitude for each BPM and then multiplying by 
the calibration scale, gave us the predicted resolutions also shown in Table [HI 

The measured resolutions for BPMs 1-2 and 9-11 are in reasonable agreement 
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Time /ms Frequency /Hz 

Fig. 14. a) The displacement of BPM 4x and b) the power spectrum of the vibra- 
tional motion as recorded by the interferometer running at 1 kHz. 

with those predicted indicating that the performance of these BPMs is limited 
by the combined electronic and digitisation errors. The resolution for BPMs 

3- 5 is significantly worse than what is predicted but is in reasonable agreement 
with the amount of vibration recorded by the interferometer (see Table |5]) thus 
indicating that the resolution of these BPMs is limited by mechanical motion. 
In order to remove this effect, we included the interferometer data in the ma- 
trix analysis as additional variables in Equation [161 This addition improved 
the resolution measurement from 0.53 to 0.45 /im (taking into account the 
geometric factor), which is equivalent to the removal of 0.35 /zm of mechanical 
vibration. This is only half of the non-rigid body motion recorded by the inter- 
ferometer. High rate (1 kHz) interferometer data indicated that the horizontal 
vibrational power spectrum peaked between 20 — 40 Hz (see Figure [H]). The 
latency between interferometer and BPM data arrival acquisition time is of 
the same order so not all of the mechanical motion could be removed. We plan 
to remedy this deficiency in future runs. 

4- 2 Calibration stability 

We took several corrector and mover calibration scans over the course of the 
dedicated 18 hours of operation in order to study the stability of the calibration 
coefficients. As the parameters of the system vary over time, these coefficients 
maintain their validity only for a limited period, after which a recalibration 
is necessary. Understanding of these effects is important for long term stable 
operation of the BPM system in the spectrometer. 

The variation of the calibration constants over the 18 hour running period 
is shown in Figures [15] and [161 The IQ phase variation is small in all the 
BPMs except for Ix and llx. The large phase change in BPM Ix was due to 
a small change in the digitiser trigger time relative to the beam arrival time. 
Other BPMs are less less sensitive to trigger fluctuations. As mentioned in 
Section [231 however. BPM Ix had significant cross-coupling from the vertical 
mode which resulted in the measured phase being time-dependent. This sen- 
sitivity produced a large change in IQ phase during the first calibration run 
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Fig. 15. Variation of IQ rotation observed over 18 hours of operation. 

and so this cahbration for BPM Ix was not used. The large changes in both 
IQ phase and scale for BPM llx are the result of the perturbation already 
discussed in Section 13. 1[ Though a new frequency was used for this section of 
the run, there is still some residual change. Consequently, we used a separate 
calibration (frequencies, IQ phase and scale) for this BPM from the time of 
the perturbation onwards. 

Apart from those exceptions noted above, the IQ phase for most BPMs in- 
dicated a drift of about 20 mrad with no significant difference between the 
mover and corrector calibration. The scale variation for the BPMs calibrated 
using the correctors was large despite the corrections of Equation [191 In the 
X direction, these variations were consistent with the statistical error (about 
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Fig. 16. Variation of scale observed over 18 hours of operation. 

±4%). In the y direction the variation was larger than could be accounted for 
from beam jitter alone. The scale variations were anti-correlated for the first 
and third stations which suggested there was a angular variation during the 
corrector scans with the pivot point somewhere between these BPM stations. 
Since only BPM 4 was equipped with a mover system, it was only possible to 
remove offset drifts, not angle drifts, using Equation [191 To measure the slope 
drift during a calibration run, we used the measured y positions of BPMs 9-11 
and BPMs 1 and 2 to calculate the overall slope of the beam and found changes 
of up to ~ 3 /irad. Over the ~ 20 m of beam line between the central BPM 
station (3-5) and the outer BPM stations, this leads to a change of ~ 60 /im 
during a cahbration step, introducing a change in scale of ~ 10 — 15%. This 
is in good agreement with the y scale variation observed. 
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Fig. 17. The temperature at the electronics racks in the experimental hall during 
the 18 hour run. 

The (9(1%) scale variations observed in the mover calibrations appear to be 
correlated to the temperature of the electronics racks (see Figure [T7|) . To 
further test whether environmental effects were behind the gain variations in 
the electronics, we applied a constant CW tone to the electronics for both the 
dipole and reference cavities and tracked drifts in relative amplitudes over the 
course of several hours. A variation of similar magnitude was found from these 
tests. 

As the variation in IQ phase was small and the large scale changes seen 
in the corrector calibrations seem to be caused by the beam drifts rather 
than electronics drifts, we averaged the coefficients obtained from the separate 
calibrations and used the mean values for the entire 18 hour run. The only 
exceptions to this were BPM Ix in which the first calibration was removed from 
the average and BPM llx where two calibrations were used, one computed 
from calibrations before the mechanical perturbation and the other from after 
it. 



4-3 Stability of BPM offsets and resolution 

We investigated the stability of the BPM system over both short and long peri- 
ods of operation. In both cases the calibration was refined using the SVD over 
the first 1000-event block of the data. The 1 hour data shown in Figure [18] is a 
zoom into the last hour of the 18 hour run (Figure [T9|) after recomputing the 
SVD coefficients. We had to select the periods of stable data-taking excluding 
the data taken during the machine tuning, big energy jumps due to klystron 
failure etc. as well as the calibration data. In addition, due to the perturbation 
of BPM llx, the SVD coefficients for BPMs 9-11 were recomputed. 

A 100 nm stability of BPMs 9-11, 500 nm of BPMs 3-5 and micron stability 
of the full system are observed over one hour of operation. The limiting factor 
for the system's stability is the stability of BPMs 1 and 2, the electronics for 
which was exposed to large temperature variations. 
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Fig. 18. The RMS and mean residual measured in different BPM stations and for 
the whole system over 1 hour of operation. Each point corresponds to 1000 events. 



For all the BPM stations, the RMS of the residual does not show any signif- 
icant change with time over the course of the 18 hour run. In contrast, the 
mean residual experiences large variations, especially for BPMs 1 and 2, and 
therefore for the linked system as well. We considered the most likely source 
of these drifts to be gain variations in the electronics. This is supported by 
Figure [2D] where we show the BPM 4 mean residual versus electronics tem- 
perature for events in which the beam position in BPM 4 was within ±25 /xm 
from the mean position over the course of the 18 hour run. Though there is a 
change in behaviour at temperatures above 27°C, below this value a significant 
correlation is seen. 
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Fig. 19. The RMS and mean residual measured in different BPM stations and for 
the whole system over 18 hours of operation. Each point corresponds to 1000 events. 
The dotted line indicates where the SVD coefficients were recomputed due to the 
perturbation of BPM llx. 



To estimate the behaviour of a triplet of BPMs when subject to changes in 
scale, we simulated three BPMs with offsets typical in our system and a beam 
orbit experiencing typical beam jitter and drifts. A scale simulated scale drift 
of 1% was implemented, similar in size to the effects seen in the experiment. 
We looked at two scenarios: only the central BPM's scale drifting only (see 
Figure [2Th) and all three BPM scales drifting by the same amount and in 
the same direction (see Figure [2Tb). In the first case, both the RMS value 
and the mean of the apparent beam residual increased gradually, while in the 
second simulation the RMS value remained almost the same, while the offset 
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Fig. 20. Residual of BPM 4x against electronics temperature. Events within ±25 fxm 
of the mean position over 18 hours of operation are shown. 




Time /A.U. Time /A.U. 

Fig. 21. Mean residuals from a simulated BPM triplet experiencing a 1% scale 
change in a) the central BPM, b) all three BPMs. 

increased gradually. The second case is consistent with the effects observed 
in our system: the physical offsets of the BPMs with respect to each other 
reconstructed by temperature-dependent gains, resulting in observed changes 
in the reconstructed offsets. 

We tried improving our stability studies for the central BPM station includ- 
ing the interferometer data into SVD computations. Unfortunately this effort 
failed as the drifts over long periods of time observed in the interferometer 
seem to be caused by the thermal expansion of the supporting aluminium ta- 
ble and the change of the refractive index of the air rather than the actual 
mechanical motion. 

We noticed that the residual shows some correlation with the position for low 
amplitude signals when the beam is close to the cavity centre (see Figure [22]) . 
which can be caused by contributions from other cavity modes. This effect is 
small, but may contribute to the observed main residual variation. 

The drifts seen in the linked system seemed to be dominated by the scale 
changes. However, over the long baseline, it is possible that variations in energy 
combined with the Earth's Magnetic Field produce a similar effect. To check 
this the energy of the beam was changed from —150 MeV to +250 MeV in five 
steps. This scan was tracked by the system as a change of the linked system 
mean residual (Figure [25]) . 
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BPM 4x Position /)j.m BPM lOx Position /^m 



Fig. 22. a) The residual between the BPM 4x position measured and predicted with 
BPMs 3 and 5 over the course of 1 hour against BPM 4x position, b) The residual 
between the BPM lOx position measured and predicted from BPMs 9 and 11 over 
the course of 1 hour against BPM lOx position. 




Relative Energy /IVIeV Relative Energy /MeV 



Fig. 23. The residual distribution of the linked system at various beam energies. 

Using a set of flux gate magnetometers, we measured the stray magnetic fields 
present around the beam line. Using the measured fields, a change in position 
of 3.1 yum in x and 0.1 fim in y is predicted for the 400 MeV energy change. We 
measured a change of ~ 3.5 fim in x with no significant change in y observed 
(see Figure [23]) which is in good agreement with the prediction. The energy 
variation over the 18 hour run was ~ ±100 MeV and so the drifts seen over 
this time period were not consistent with a change in energy. 



5 Conclusions 



We have successfully commissioned eight cavity BPMs of differing designs and 
properties divided into three BPM stations in the End Station A beamline 
at SLAC. The first BPM station consisted of two rectangular cavity BPMs 
originally designed for use in the A-line. They demonstrated resolutions in 
X and y of 1.1 fim and 2.2 fim respectively and were stable to ±5 fim over 
one hour while drifting by ±10 /im over 18 hours of operation. The resolution 
was limited by digitiser errors and noise. The drifts can be explained by the 
changes of electronics gains. 

The second BPM station consisted of three prototype ILC linac cylindrical 
BPMs with the central BPM mounted on a dual-axis mover. These demon- 
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strated resolutions of 0.53 /xm in x and 0.46 iim. in y. Mechanical motion was 
found to dominate the resolution. These BPMs were stable to ±0.25 /im over 
a period of 1 hour and ±1 fim. over a period of 18 hours. The stability was 
influenced by low amplitude effects and mechanical vibration on short time 
scales and scale changes of ~ 1% magnitude over long time scales. 



The final BPM station in the beamline consisted of three rectangular cavity 
BPMs, originally designed for use in the SLAC hnac. They demonstrated 
resolutions of 0.19 /im in x and 0.17 /xm in y with a stability of ±50 nm over 
the one hour period. Long term stability of this station could only be measured 
over 10 hours due to a perturbation that altered the calibration of BPM llx. 
However, over these 10 hours, the triplet achieved a stabihty of ±500 nm in 
X and ±750 nm in y. As with BPMs 3-5, we think the stability was limited 
by low amplitude effects and scale drifts, but there was no information on the 
mechanical motion available for this triplet. 



When combining all the BPM stations to measure the precision of the orbit re- 
construction over the whole baseline, a resolution of 0.82 /im in x and 1.19 /im 
in y was achieved. The system was stable at the micron level over the course 
of one hour. The long term stability was affected by relative scale drifts across 
all the BPMs therefore drifts of the order of ±10 jjim were observed over 18 
hours of operation. 



In order to improve the system we had in 2006 and be able to perform full 
tests of a magnetic spectrometer prototype, we have added several upgrades 
to the beamline and BPM systems. Four steel core dipole magnets have been 
installed to form the magnetic chicane with BPMs 3 and 5 now measuring 
the incoming beam position, BPMs 9-11 measuring the outgoing beam po- 
sition and BPM 4 having been moved to measure the beam position at the 
mid-chicane location. To improve the stability of the BPMs, we have added 
a periodic sine wave calibration tone. This is applied through the BPM elec- 
tronics at a rate of ~ 0.1 Hz and allows continuous gain and phase monitoring. 
This should allow us to track calibration changes online and correct for any 
variation. The interferometer has been extended to measure mechanical drifts 
between the new incoming BPM station (BPMs 3 and 5) and the mid-chicane 
BPM. Hclmholtz coils have been installed to quickly dither the beam and 
therefore reduce the dependence of the calibration procedure on beam drifts. 
We have commissioned these upgrades during 2007 and are planning to con- 
tinue in 2008, so we will report on the results in our next publication. 
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